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The technique of quantum repeaters is a promising candidate for sending quantum states over long 
distances through a lossy channel. The usual discussions of this technique deals with only a finite 
dimensional Hilbert space. However the qubits with which one implements this procedure will "ride" 
on continuous degrees of freedom of the carrier particles. Here we analyze the action of quantum 
repeaters using a model based on pulsed parametric down conversion entanglement swapping. Our 
model contains some basic traits of a real experiment. We show that the state created, after the 
use of any number of parametric down converters in a series of entanglement swappings, is always 
an entangled (actually distillable) state, although of a different form than the one that is usually 
assumed. Furthermore, the output state always violates a Bell inequality. 



I. INTRODUCTION 

Entanglement cannot be created by local operations 
and classical communication between the parties. How- 
ever it was shown in Ref . (see also @) that there exists 
an operational scheme, such that particles can get entan- 
gled without ever having interacted in the past. One of 
the intruiging things about this phenomenon, which has 
been called entanglement swapping, is that it shows that 
one cannot always tell whether particles are entangled 
by looking at their "common history" . Or the concept of 
"common history" must be suitably enlarged. Note that 
entanglement swapping process is a specific case of quan- 
tum teleportation |3|. The first experimental realization 
of entanglement swapping was reported in 0. The ex- 
periment was a direct realization of the experimental pro- 
cedure given in 0, and modified in 0. 

As a simple illustration of this phenomenon, consider 
the situation in which Alice and Bob, share the singlet 
= -^7|(|01) — 1 1 0) ) , and Alice shares another singlet 
with Claire. Alice now makes a projection measurement 
on her parts of the two singlets in the Bell basis, given 
by the states 



|^=) = ^(|00)±|11», 
h/-±> = ^(|01>±|10». 



(1) 



It is easy to check that if Alice now communicates (over a 
classical channel) the result of her measurement to Bob 
and Claire, they will know that they share one of the 
Bell states given by eq. Q). Note that in principle the 
particles of Bob and Claire may not have interacted in 
the past, although they share entanglement after Alice's 
classical communication to them. 

Apart from this fundamental perspective, entangle- 
ment swapping is also important in quantum communi- 
cation applications. When sending a quantum state over 
a noisy channel, the probability that it reaches the recip- 
ient, decreases with the length of the channel. However, 
the detectors at the recipient's end are usually (rather 
invariably) noisy, and this noise is independent of the 
length of the channel. Thus after a critical length, the 
signal is useless. To circumvent this problem, a proposal 



was provided [|| that places a number of nodes in be- 
tween the sender and the recipient of the signal. Entan- 
gled states are first shared in these shorter segments (i.e. 
between all successive nodes) and thereafter distilled 
to obtain highly entangled states between all successive 
nodes (Fig. Finally entanglement swapping is carried 
out at all nodes to obtain highly entangled states between 
the sender and the recipient (henceforth called Alice and 
Bob respectively). It was shown |(j that this procedure, 
called quantum repeaters, would lead to highly entangled 
states between the ends of a noisy channel of arbitrary 
length, with only a polynomial increase in time and log- 
arithmic increase in local resources. 

Usual discussions on quantum repeaters deal with only 
a finite dimensional Hilbert space. But the qubits with 
which one implements this procedure will "ride" on con- 
tinuous degrees of freedom of the carrier particles. In 
this paper, we address the problem of implementation of 
this procedure of quantum repeaters, for the entangled 
states prepared between successive nodes by spontaneous 
pulsed parametric down conversion . The discussion 
of the process will follow the ones of Refs. 00. The 
experimental realization of entanglement swapping fully 
confirmed the validity of this description . Actually in 
the experiment, polarization entanglement was utilized. 
But it is elementary to show the equivalence of such ex- 
periment with ones involving path entanglement, which 
will be our model here (see Ref. |8|). In this paper, we 
assume that the noise in the channel from a parametric 
down conversion crystal to the nearest nodes is negligible. 
Entanglement swapping is carried out at all the nodes. 

The description of entanglement swapping that we con- 
sider in this paper, is still a toy model. Nevertheless it 
contains some basic traits of a possible real experiment. 
We show that the final state prepared between Alice and 
Bob is a so-called maximally correlated state, which is 
always entangled (actually distillable) and always vio- 
lates a Bell inequality. For a wide range of pulses and 
filters, including Gaussian pulses and filters, the output 
two qubit state created between Alice and Bob turns out 
to be a mixture of two Bell states, 
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FIG. 1: Quantum repeaters: A chain of maximally entangled states obtained after distilations, are shared between the sender, 
the nodes at A, B, . . ., N and the recipient. Bell measurements are performed at the nodes. This results in an entangled state 
shared by the sender and the reciepent. 



where j^) are given by eq. (QJ and where < V < 1. 



II. DOUBLE ENTANGLEMENT SWAPPING 

Consider the case of three parametric down conver- 
sion crystals producing three entangled states between 
the sender (Alice) and A\, A2 and Bi, and B2 and the 
recipient (Bob) (see Fig. |2J0]. 
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FIG. 2: Quantum repeaters: A chain of three states. 

Entanglement swapping will be carried out at A1A2 
and B\B2- In Fig. we give a schematic description of 
entanglement swapping in an array of three type I spon- 
taneous parametric down conversions (PDC) [l], |5| • Note 
that frequency filters are in front of every beamsplitter, 



in the "internal" part of the device. They are necessary 
to make the photons emerging independently from two 
different sources, indistinguishable by the detectors be- 
hind the beamsplitters (for the physical reasons for this, 
see 0,13). 

We make the simplifying assumptions that the optical 
lengths of all source-detector paths are equal and that 
phase shifters work in the range of the order of the wave- 
length (i.e. between and 2tt). This enables us to neglect 
all retardation effects. 

The description of the two-photon initial entangled 
state will depend on whether the corresponding PDC is 
an "external" or an "internal" one, in the series of PDCs. 
In any series of PDCs, there will be two external ones, 
while the rest will be called internal. For example, in the 
case of three PDCs, as described in Fig. |3 PDC-I and 
PDC-III are externals, while PDC-II is an internal one. 
If the "idler" photon emitted by the external PDCs in 
Fig. |3| (produced by a single pulse from a laser pump) , 
manages to pass via the filters, the resulting two-photon 
state is given by (see Appendix) 



ip ext (x,y,x ,y ) 



duji 



duj s / duj p F(uj p )A(uj p -u>i - w„)/(a;i)(<4 (ui)a\{w s ) + a*,(a; i )aj,(w s )) |fi) , (2) 



where oJ(w») (al(u) s )) and a^(wi) (a 1 ,^)) are the ere- at PDC-I is 



ation operators of photons of the idler (signal) of fre- 
quency u>i (uj s ) respectively in beams x and x {y and 
y ). The function F represents the spectral content of 
the pump pulse of frequency uj p and / is the transmis- 
sion function of the filters and is assumed to be centered 
at u)q/2, where in turn, luq is the central frequency of 



the laser pump pulse. The function A(a> p 



Uj) is 



due to the phase-matching condition of the PDC pro- 
cess. | il) is the vacuum state. The state produced 



J 



tpext(b, a,b , a )) and that at PDC-III is 

ipext{fi e if i e )^- We assume, in all our considerations, 
the perfect case and so we will replace our A by the Dirac 
delta function. Here and henceforth, unless stated oth- 
erwise, we ignore normalization of states. 

If the two photons produced by the internal PDC 
(PDC-II) in Fig. [21 manage to pass the filters, their state 
acquires the following form (see Appendix): 



ipint(d,c,d ,c)\ = / dJi I du)' s I duj' p F((j' p )A(uj' p - J i - wl)/(a;-)/(a; s )(4(w-)4(w s ) + a^(w-)a^, (w s )) |Q) . (3) 
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FIG. 3: Entanglement swapping in a chain of three parametric down converters. The boxes denoted as PDC-I, PDC-II, and 
PDC-III symbolizes the full parametric down conversion device, including the pump laser, crystal, a set of apertures that select 
suitable output directions, as well as suitably oriented mirrors that direct the radiation into the interferometers shown in the 
figure. The radiation at the detectors ii, ...,14 and i 4 is observed by detectors. The interferometers are built out of 

mirrors M and 50-50 beamsplitters (BS). We need to put suitable frequency filters in front of all BS leading to the internal 

outputs and y The symbols 4> an d 4> denote the local phase shifters. For simplicity, we shall study the case when only 
four internal detectors (the "i"s) fire (those marked with asterix). Due to the specific properties of downconversion 0,13, no 
filters are needed in the output beams of the device, a and /. 



Note that there is an extra filter function in this case. 
This is a signature of the fact that both the photons from 
PDC-II are detected by the internal detectors and to 
reach there, they must pass the filters. In the general case 
of an arbitrary number (say, n) of PDCs in a series, there 
will be two entangled pairs described by eq. J2J), while 
there will be n— 2 "internal" entangled pairs described by 
eq. Note here that previous works on entanglement 
swapping with two PDCs |l|,|5j did not need to consider- 
any "internal" entangled pairs. 

In the case of entanglement swapping by 3 PDCs (as 
shown in Fig. |2J), our initial state, if all photons manage 
to pass the filters, is given by 



i^extU \e,f ',e')\ ip int (d,c,d \c)j ip ext (b,a,b' ,a')^> 



Now suppose that the detector i\ fires at time t\, i 2 at 
time t2, «3 at time £3, and i± at time t^. Then the wave 
packet collapses into the state 

I*! h, h, t 3 , U) = il(*l)t 2 (t2)*3(*3)*4(t4) I* ) (4) 



where, for example, 

ii{t\) = J duj exp{— iujti}ai 1 (w). 

Note that the scalar product of J dio &L\>\iujt\\a\^ (u) 
with a single photon state 



gives the probability amplitude to detect this photons at 
time t\. Let us assume that our 50-50 beamsplitters (BS) 
are the symmetric ones. That is, one has, e.g. 



a\ M = ^ [ a t> M + ia l> ( w } 



Note that the creation operator of the reflected beam 
always enters the relation for aj with an i factor. As 
we consider here the idealized case, we substitute A by 
Dirac delta functions to obtain 



I*;ti,t 2 ,t3,t4) - J dtF*(t)f{t-U)f{t-ti)a}[g,t 3 ]al[g,t 2 }+ J dtF*{t)f{t-U)j{t~t 2 )aJ fl [gMW al [gM], (5) 

I 

where we denote for example, with 

g(u a ,t 3 ) = J dtexp{tuj s t}F*{t)f{t-t 3 ). 
a f [9>h] = / duj s g(uj s ,t 3 )a f (uj s ) After the photons in the state |^ ; ^, t 2 , t 3 , t 4 ) pass 



4 



through the phase shifters (as depicted in Fig. |3J, the 



new state, 'J'; t\, t 2 , t 3 , t 4 , (f>, 4> y, is obtained by replac- 
ing a\ by exp{i0}a^, and at by exp{i(j> }a|. Let us as- 
sume now that the photons (after passing through phase 
shifters and 50-50 beam splitters (cf. Fig. |3J|) are de- 
tected at a + and / + at times t a and tf respectively. 
Hence the amplitude of such a process is 

A a+ J + {t a , tf, tl, t 2 , £l) = 



(/+(*/), a+(t a ) I 



^',tl,t 2 ,t 3 ,t 4 , 



where, for example, 



\f+(tf),a+(t a )) = f dwexp{wtf}a f+ (ui) 

J dco cxpjicj t a }a\ + (cu ) |fi) . 

Writing it out explicitly, one obtains 

Aa+,f+ = exp{i(0 - )}Ti + T 2 , 

where 

?l = T(t Q , tf,ti,t2,ts, t 4 ) 
= J dtdt dt J duj s duj s exp{i(o; s t + uj s t )}x 
exp{-i{ujj f + u s t a )}F* {t)F* (t )F* (t" ) x 
f(t-t 3 )f(t -t 2 )f(t" -U)f(t" -ti) 



(6) 



and 



?2 — r^o, t/, f 2 , tl, t4j t3)- 



However, due to the finite time resolution of the detec- 
tors, the precise detection times are not known. There- 
fore the full probability, P, of the process is obtained 
after integration of the square of the amplitude over the 
detection times: 



P(a+,f+) 
+2 cos((A - 



J[\Ti 



) ) I T\ T£ I ] dti dt 2 dis dt^dtf dt a 



(7) 



(where cj>o is a phase). We integrate here from —00 to 
+00, because the time resolution of detectors is usually 
by orders of magnitude bigger than the duration of the 
pumping pulse. We can now calculate the interferometric 
contrast, or the visibility 



V* = 



max P — min P 
max P + min P ' 



of the two particle interference process as observed by 
changing the values of the local phase shifts in the two 
external interferometers. (The index 3 stands for the 
number of PDCs involved.) The resulting value is given 
by 



V* = 



2 J \T 1 T^\dt 1 dt 2 dt 3 dt 4 dt f dt a 

J(\Tl\ 2 + l^l^dhdhdtsdtidtfdta 



(8) 



We will now write down the explicit forms of the ex- 
pressions in the numerator and denominator of V3. We 
have 

-F1234 = / \T1T2\dt1dt2dt3dt4dtfdta_ 
= J dt^dt^dt^dt^ J dtdt dt dt" x 
\F(t)\ 2 \F(t')\ 2 F*(t")F(F)x 
f(t - t 3 )f*(t - U)f(t - ta)/*(t' -Ji)* 
f(t" - t 4 )f(t" - t x )f*{t" - t 3 )f*(t" - t 2 )|. 

(9) 

Throughout the paper, a "bar" will represent a new vari- 
able. This should not be confused with complex conju- 
gation, which we denote here by a "*". 
A similar simplifying gives 

G1234 = J \Ti\ 2 dtidt2dt 3 dtidtfdt a 

= J dtdt dt dt" J dtidt 2 dt 3 dt4X 
\F{t)\ 2 \F{t)\ 2 F*{t)F{T)x 
\f(t-t 3 )\ 2 \f(t_-t 2 )\ 2 f(t" -t 4 )f(t" 
/*(*"-**)/*(*" 



Note that in G1234, the indices are ordered, 
easily find that 

J \T 2 \ 2 dtidt 2 dt 3 dt 4 dt fdt a = G 2 i i3 . 



(10) 
One can 



(11) 



Hence the visibility of the two particle interference that 
can be obtained due to the two-fold entanglement swap- 
ping (Fig. |2J), is given by 



V„ = 2Fl234 

G1234+G21 



(12) 



III. ARBITRARY NUMBER OF SWAPPINGS 

Consider now a chain (i.e., arranged in a series) of 
any number, n, of PDCs (see Fig. Let aW and 

a W denote the paths within the external interferometer 
into which one of the photons from the first PDC en- 
ters. The paths within the other external interferometer, 
into which enters the photon from the last PDcs will be 
denoted as and a W. Let wi 1 -* and lo^ be the cor- 
responding frequencies of the photons. If all the photons 
manage through the frequency filters, the state is 



(1) 







(13) 



with the states li/' ) an d IV' ) of the type character- 
istic for the external PDCs and (I = 2, . . . , n - 2) 
for the internal ones (compare eqs. and J3J)). Imagine 
now that again all internal detectors ik, k — 1,2, ... ,2n— 
2 fire. The final state into which the state in eq. i|13|) 
collapses is of the following form: 

\ Xt ) = J^^" ) (X t (^ 1 \^)at ( „ ) (^Vl' (1) (^ 1) ) 



)at, ( „ ) (4" ) )ai( 1 )(^ 1) ))l«) 



(14) 
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The functions X t and Y t depend on the detection times at 
the internal detectors (along with the frequencies). The 
subscript t stands for the full set of times of detection at 
the internal detectors, that is t = ti, t%, . . . , t2n-2- To see 
this structure of the state for n = 4, see Fig. 0] 

The effective mixed state, p swa P, that we obtain, is an 
incoherent sum of the state Xt , with the sum (integration) 
being over the detection times: 



dt\ X t)(Xt\- 



where dt = n^dif. 

Again we now assume that optical lengths of all source- 
detector paths are equal and that phase shifts <fi an d <t> 
are of the order of the wavelength (i.e. between and 2tt). 
Since the two output photons are fully distinguishable 
(i.e. their sources are known), it would be no harm to 
abandon the second quantized description. Thus, one can 
rewrite the formula (|14|l in the following way. Note that 
{uS) |f2) is a single photon state. The photon is in 
path and has frequency lu. Therefore from the first 
quantized point of view, this state can be replaced by 
the tensor product \a) \lo), where \a) describes the path 
variable and \oj) the frequency variable (energy). Using 
such notation, eq. I|14|) can be put as 



and similarly for 



,(„) 



\ Xt ) = JdU 1} dJ s n) x 



,(«) 



.(«) 



(15) 

The frequency degrees of freedom in p swa P can be traced 
out to obtain the (unnormalized) state of the path (i.e. 
the qubit) degrees of freedom as 



where 



(IX^iV^)! 2 a' (1) )(a' (1) ®|oW>(oW 
+\Y(U 1 \ui n) )\ 2 |a (1) ) <a (1) | <8> |a' (n) \ 
+ [X(J s 1 \U n) )Y*(U 1 \U n) )x 
a'Wj>(a«| <g> |a(™)) (a'W +h.c.]). 



Here h.c. denotes hermitian conjugate of the term before 
it within square brackets. 

After normalization, one gets 



-CT = w-d b loo> (00| + c |11) (11| 

+ [a|00)(ll| + h.c.]) 



(16) 



Here a = J dwdtXY* , b = J dudt\X\ 2 , c = J dujdt\Y\ 2 , 
and 



|0) = |a«), |l) = |a'( 1 )) 



) and 



Writing a 



r a exp{i9 a } (r a > 0, 6*0 real), and redefining the state 



|1) of the first particle as exp{— i9 a } 
ized state p™?? reads 



W ), the normal- 



path 

swap 
Ppath 



^(6|00) (00| + C |11) 
[r„|00)(ll|+h.c.]) 



(HI 



where r a , 6, and c are all positive. 
One can now check that 



Vn 



2r a 
b + c 



(17) 



(18) 



where V n is the visibility of the two-photon interference 
in the external interferometers (like those in Fig. 

Let us give the values of the parameter in eq. (|18|l for 
the case of Fig. [3] One has 

r a = J du) s duj s J dtidt2dt3dt4 J dtdt dt J dtdt'dt"x 
exp{i(ujj + u s t')}F* (t)F* (t )F* (t ) x 
/(* - t 3 )f(t -h)f(t" - U)f_{t -h)x 
exp{i(u'^t +_u s t')}F%t)F* (t' )F* (t" ) x 

f(t - u)f(7 - tt)f(t" - t 3 )f(t" - 1 2 ). 

Comparing with eq. 10 , one can verify that 
r a = F 123 4- 

Similarly using eqs. (|10|l and 1)11(1 . one finds that b and 
c are respectively G1234 and G2143. Using eq. (|12[) . one 
therefore obtains the relation in eq. I|18l) . for the case of 
three parametric down conversions. It is straightforward 
to see that the same relation holds for an arbitrary num- 
ber of PDCs. As we have mentioned earlier, the case of 
two PDCs is slightly different, in the sense that there are 
no internal PDCs. However it is easy to check that the 
relation in eq. ((18JI is true even for the case of two PDCs. 

Whenever r Q ^ 0, the state /0*^f has nonpositive par- 
tial transpose 03 • Therefore the state is entangled when- 
ever r a ^ ^l| (see also 0]). The entanglement of 
formation of x P ath is El 



£ = H(±(l + y/T=V*)) 



(19) 



where H{x) — — aTog 2 x — (1 — x) log 2 (l — x) is the bi- 
nary entropy function. Note that the visibility V n is the 
so-called concurrence 0| of the state p^th '■ ^ n ^ act ' 
entanglement of formation is known to be additive for 
the state p™y? (see also E3)> an d hence the ex- 

pression displayed in (|T§|> is also the entanglement cost 
(asymptotic entanglement of formation) of pt at ^ '■ States 
with nonpositive partial transpose in 2 £g> 2 are distillable 
\vl\. Therefore the state pV^^, being in 2 ® 2, is also 
distillable whenever r a 7^ 0. 
Bipartite States of the form 
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FIG. 4: An exemplary extended chain of swappings. Imagine that the initial states emitted by the four sources were entangled 
two photon states, with the property that the one photon is emitted to the right and other one to the left, and the emissions 
are such that both photons are in primed beams or both are in the unprimed ones. Now if the detectors ii,...,ie fired (the 
ones with asterix), one has only two consistent possibilities, (i) i\ fired due to source 1, but then in fired due to source 2, 
therefore 13 fired again due to source 2 whereas ii fired thanks to source 3, thus is fired also because of source 3, and finally 
ie fired because of source 4. This implies that a photon is in the beam a' 4 ' and in the beam a (ii) The second possibility 
is that the ii firing is due to source 2 and etc. Thus the other element of the final superposition must consist of photons in 
beams o and in the beam a' 4 '. 



are called maximally correlated states [l8j. It is known 
that for such states in 2 eg) 2, the entanglement cost is 
strictly higher than its distillable entanglement [lH Il5| 
(cf . |18l Il9| ) . We therefore have irreversibility in asymp- 
totic local manipulations of entanglement for such states. 
The state pp^f being a maximally correlated state in 
2 (g) 2, would also have its entanglement cost strictly 
greater than its distillable entanglement. 

Additionally, the state p^jf violates a Bell inequality 
20 1 whenever r a 7^ 0. The maximal amount of violation 
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2 
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Vn 



(20) 



The r a = case corresponds to null visibility and can 
therefore be ignored. We plot the entanglement of for- 
mation (which in our case is entanglement cost) and the 
maximal amount of violation of Bell inequality against 
the visibility V n in Fig. [3] We therefore obtain the im- 
portant fact that the output state, resulting from a se- 
ries of entanglement swappings, is always entangled and 
always violates local realism. Earlier works on entangle- 
ment swapping process via two parametric down convert- 
ers pi made the tacit assumption that the output is 
a pure state admixed with white noise, and reached the 
conclusion that the swapped state is separable and does 
not violate local realism for low visibilities. 

We will now discuss an interesting point with respect 
to violation of Bell inequalities of the state pj^f ■ Let 
us, for this state, find the T-matrix, whose elements for 
any state p, of two qubits, is given by 

Tij = tv{o-i ® <TjQ), (i,j = x,y,z). 



In the above formula, we treat our qubits formally as 
spin-1/2 particles. For the state Pp„tT> Txx — Vn, T yy — 
—V n , T zz = 1 and the rest are vanishing. A necessary 



FIG. 5: Plot of entanglement cost (£), given by eq. I19L 
and the maximal amount of violation of Bell inequality (B) , 
given by eq. I|2()|l . against the visibility V n for the state Xpath 
obtained after the detection of the idlers in the entanglement 
swapping process with n parametric down conversions. 



and sufficient condition for a state of two qubits to violate 
local realism, for two settings per site, in the plane of n 
and n is given by j2^ 



E 



TA > 1. 



In our case (i.e. for the state p™th)> T ij = 2V n, 

so that for V n > -7=, one obtains a violation of local 

V 2 

realism in the y plane. In the case of qubits defined by 
the output paths of the multiple entanglement swapping 
devices considered here, x — y plane spin observables are 
equivalent to the measurement in the output ports of the 
external interferometers. That is, the n-a operator for an 
n in the x — y plane, is equivalent to a device consisting 
of a phase shifter in front of a 50-50 beamsplitter and 
two detectors behind it. 

The x — y plane does not provide a violation of local 
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realism for < V n < although the state is still entan- 
gled (and distillable) in that region. For these lower val- 
ues of the visibility V n , we have to consider other planes 
for obtaining a violation. For example in the x — z plane, 
52i4=*,* T ij = 1 + V n> fOT the state p s p ^, and thus a vi- 
olation is always obtained in this plane. Hence, for these 
lower visibilities, one must use other measurement planes 
to obtain a violation. Such violation can therefore be ob- 
tained, only by the Mach-Zehnder interferometers at the 
output ports of the entanglement swapping device (Fig. 
HjJ). This is due to the fact that such a device is capable 
of performing any U(2) unitary transformation. 



One has to bear in mind that the above results were 
obtained under the asuumption that we deal with perfect 
entanglement swappings, especially no noise was allowed. 
It is obvious that for sufficiently low V n (see Fig. |SJ, both 
entanglement and violation of Bell inequalities would dis- 
appear, even for very minor noise admixtures. 
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FIG. 6: The schematic diagram shows that to obtain better 
violation of local realism, the output ports of the entangle- 
ment swapping (contained in the square box in the figure) 
should be put to the Mach-Zehnder interferometers. 

For a wide range of pumps and filters used in the swap- 
ping process, one will have 



For example, this is the case when the pulse and filter 
functions are Gaussian, i.e. 



and 



77, / v r (Wp-^o) 2 

F(u p ) = exp{ p — } 



tt \ r (^p-^o/2) 2 -, 
/(w p ) = exp{ 1—^ }. 



2(7, 



(21) 



(22) 



In such cases, the normalized state created after the 
idlers have been detected in the multiple entanglement 
swapping process with n PDCs is 

PTatn = |(|00)(00| + |11)(11| 

+ K (6 - c) |00)(ll| + K (b - c) |ll)(00|), 

which can be rewritten as 



swap 

Ppath ~ 



1+V, 



(b=c) 



\<t> + )(ct> + \ 



1 - V, 



(6=c) 



2 / \~ 2 

where j^) are the Bell states given by eq. 



4>-)(<t>- 



APPENDIX A: THE TWO-PHOTON STATE 
PRODUCED BY PDC 



This appendix is essentially meant to provide a 
"derivation" of the state (displayed in eqs. J5J and (J3J) 
produced in a parametric down conversion. It is a re- 
formulation of the theory of parametric down conversion 
which can be found in, e.g., pjj.!^. 

The phenomenon of spontaneous parametric down 
conversion (PDC) is a spontaneous fission of quasi 
monochromatic laser photons into correlated pairs of 
lower energy. All that takes place within a non-linear 
crystal. The probability of a single laser photon to fission 
is very low, but in a strong laser beam, the frequency of 
the phenomenon is quite high. The new photons, custom- 
arily called "signal" and "idler" have some basic proper- 
ties. First of all, the wave vector ko of the laser photons 
is related to those (fcj and k s ) of the idler and the signal 
by ko w ki + k s . One has to stress that this relation holds 
within the crystal (and can be thought of as an approx- 
imate conservation of the linear momentum) . Secondly, 
the frequencies ojq, oji, and cu s of the laser photon, idler 
and signal satisfy oj ~ u>i + lu s . And finally, the emerg- 
ing pairs of photons are highly time-correlated. That is, 
if their optical paths from the source to the detectors 
are equal (which we assume in this paper), the detection 
times are equal too (up to the time resolution of the de- 
tection system). The relations between the wave vectors 
and between the frequencies are called phase matching 
conditions. 



1. Crystal-field interaction. 

Let us recall that in the interaction Hamiltonian of 
the electromagnetic field with an atom or molecule, the 
dominating part is H a -f ~ jl e ■ E(x,t). That is, it is 
proportional to the scalar product of the dipolc moment 
of the atom or molecule with the local electromagnetic 
field. Now the electric polarization of a material medium 
is given by the mean dipole moment of the atoms or 
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molecules from which the medium is built (per unit vol- 
ume). Let p(x,t) stand for the local polarization of the 
volume 8V , which contains the point x (the volume 5V 
is macroscopically small but microscopically large). The 
principal term of the interaction Hamiltonian for crystal 
and field must have the form 

H mt ~ f p(x,t) ■ E(x,t)d 3 x, (Al) 
Jv 

where V is the volume of the crystal. From the micro- 
scopic standpoint, the above formula reads 

flint ~ f '%2iin8 ( ®@-3e n )-fi(Z,t)dPx, (A2) 
^ n 

where x n is a symbolic representation of the position of 
the n-th atom endowed with a dipole moment fl n . The 
summation is over all atoms (the word "atom" standing 
for any stable aggregate of charged particles, like atoms 
themselves, or ions, or molecules) of the medium. One 
can see that the formula l)A2|) agrees with fl a -/- If one 
introduces the averaged (macroscopic) polarization (av- 
eraged over the macroscopically small volumes SV), we 
get I|A1I) . One can assume that E(x, t) interacts with 
p(x, t) only in the point x, thus the i-th component of 
polarization is in the most general case given by 

Pi(x,t) = Ef=iXy (x)Ej(x,t) 

+ E J >ixgm(^^(£,t) + - ) 

where x\j are X«i are ^ ne (macroscopic) polarizability 
coefficients. They are in the form of tensors. This is 
due to the fact that polarizablility may depend on the 
polarization of the incoming light. Note here that for 
any crystal which is built of molecules which are centro- 
symmetric the quadratic term of the polarizability van- 
ishes. Therefore the effect exists only in birefringent me- 
dia. 

In the case of a perfect "nonlinear" crystal, we assume 

(2) 

that Xijk\%) nas the same value for all point within the 
crystal. The nonlinear term of the polarization gives the 
following term in the interaction Hamiltonian (cf. ilAlfl ): 

Hint ~ J v p(x, t) ■ E(x, t)d 3 x 
= f v P lin &t) ■ E(x, t)d 3 x 
+ J v P non (x 1 t)-E(x, t)d 3 x, 

where p hn (p non ) is the linear (nonlinear) term of polar- 
ization. The nonlinear part of the Hamiltonian is 




The quantized field can be expressed (in the interaction 
picture) as 

E(x,t) = VL., f d 3 k-^ X 

e(k,p)a~ (k,p) exp{i(k ■ x — ojt)} + h.c. (A5) 
= £(+)(.?, t) + E(~\x, t), 



where £j' \x, t) = [E^ (x, t)]\ and the summation is 
over two orthogonal linear polarizations, h.c denotes the 
hermitian conjugate of the previous term, e(k,p) is a 
unit vector defining the linear polarization. The symbol 
a~{k,p) stands for the annihilation operator of a photon 
of a wave vector k and polarization e(k,p). The principal 
commutation rule for the creation and annihilation op- 
erators is given by [a(k,p), a?{k',p')] = 8 p yS^ 3 \k — k'), 
[a^(k,p),a^(k',p')] = and [a(k,p), a(k',p')] = . As we 
are interested only in the PDC process, we will neglect 
the depletion of the laser field and assume that the total 
field within the crystal is E(x, t) = E Laser (x, t)+E q (x, t), 
where the laser beam is described by a classical electro- 
magnetic field E Laser . In reality, the laser field is a mix- 
ture of coherent states, and its phase is undefined, but 
this is of no consequence to us here. The field E q is de- 
scribed in a quantum-electrodynamical way. It describes 
the secondary photons emitted by the crystal. The down 
conversion takes place, thanks to the terms in (|A4I) of 
the form J v J2i J kX^kH[ aser E J E k d 3 x. Only those terms 
of Ej and Ek which contain the creation operators, can 
give rise to a two photon state, after acting on the vac- 
uum state 1 17). The creation operators can be found only 
in the so-called negative frequency terms of the electro- 
magnetic field operators (that is, in those which contain 
the factors exp{— i(k ■ x — ujt)} (cf. EH))- L et us there- 
fore forget about all other terms and analyze only 




For simplicity let us describe the laser field as a 
monochromatic wave EL aser (x,t) — xEo(exp{i(ko ■ x — 
Loot — </))} + c.c), where 2Eq is the field amplitude and 
c.c denotes the complex conjugate of the previous ex- 
pression. (Since an arbitrary electromagnetic field is a 
superposition of the plane waves, it is very easy to get 
the general description.) Then from (|A6p . one gets 

h nl ~J v EM^ E °x 

[exp{i(/co • x — u>ot — (f>)} + c.c] x 

E P / d3k f(^)Kk,p)aHk,p) exp{-ifc • x}x 

E P ' / d3k ' /(w)e(fc,p)a f (fc',p) exp{-ik> ■ x}d 3 x 
+h.c, 

(A7) 

with f(ui) being a factor dependent on u). Its specific 
structure is irrelevent here. Extracting only those el- 
ements of the above expression which contain x, one 
sees that their overall contribution to (|A7|) is given 
by J v d 3 x exp{ix(±/co — k — k')}, which we write as 

A(fco — k — k'). If we assume that our crystal is a cube 
Lx Lx L, then for L — > oo, A approaches the Dirac delta 
5(±fco — k~k r ). Thus, we immediately conclude that the 
emission of the pairs of photons is possible only for the 
directions for which the condition ko sa k + k' is met pBj . 
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Eq. (|A7|) can be therefore put in the form 

H NL ~ E P y Sfk J d 3 fcA(fc - * - k>)A e p f J, x 
exp{— iuJot}(v (k,p)a) (k' ,p') + h.c, 

(A8) 

where A e p f J, = Y^j,k E oxf] k ^^P)h{k' ,p') is the effec- 
tive amplitude of the laser pump field (which serves as a 
laser-crystal coupling strength). This Hamiltonian fully 
describes the basic traits of the phenomenon of down 
conversion. In the so called type II down conversion, the 
laser pump beam is ordinary wheras the idler and sig- 
nal photons are extraordinary. Thus we shall replace the 
general A e Jf, by Fq. 

1. The state of photons emitted in the PDC 
process. 

We are interested in the process of production of pairs 
of photons. Therefore we shall assume that the pump 
field is rather weak, so that the events of double pairs 
emissions are very rare. Therefore one can use the per- 
turbation theory. 

The evolution of the state | *£>(£)) (in the interac- 
tion (Dirac) picture) of the photons emitted in the 
PDC process is governed by the Schrodinger equa- 
tion ih±\* D (t)) = H% L (t)\^ D (t)), where H% L (t) = 
exp{±-H t}H NL exp{-j-H t}, and H +H NL is the total 
Hamiltonian of the system. Therefore 

\* D (t)) = \v D {t )) + ±ftm L {ti)\* D {t>))dt> 

where we have replaced the time dependent state in 
the integral by its initial form | \& £>(*(>)) using the first 
order of the perturbation calculus. In the Dirac pic- 
ture, the creation and annihilation operators depend on 
time as a D (k,p,t) = exp{iujt}af(k,p) and ao(k,p,t) = 
exp{iujt}a(k, p) . In we put to = 0, and we take the 

vacuum state (no photons) |fi) as the initial state |^(0)). 
We are interested only in the term with the integral, be- 
cause it is only there that one can find creation opera- 
tors responsible for the spontaneous emission of pairs of 
photons. The photons interact with the laser field only 
during the time of the order of — . The interaction sim- 
ply ceases when they leave the crystal. Therefore, as the 
annihilation operators when acting on the vacuum state 
give zero, one can write 

E P g> I d 3 kj d 3 k'F A(k -k- k')A(uj + uj' - uj q ) 
a^k,p)aHk',p')\n), 

L_ 

where we have written J " dt' exp{it'(uj + uj — luo)} as 
A(ci> + uj' — (jj ). As L — > oo, A behaves as a Dirac delta. 
Thus the allowed frequencies of the emissions satisfy the 
relation [2(| ujq ~ uj + uj' ■ Since — uj, —uj' and — ujq are 
typically of the order of 10 4 the function A(uj + uj' — ujq) 
is very close to 5 (to + uj' — ujq). 



3. Directions of emissions 

We know that uj = \k\ , c where , c . = c(ui.p) is 
the speed of light in the given medium, which depends 
on frequency and the polarization. Using this relation as 
well as phase matching condition for frequencies, we get 
the condition \k \c(uj ) ~ \k\c(ui,p) + \k'\c(ui' ,p'). There- 
fore, the emissions of the pairs are possible only when the 
phase matching and the above condition are both met. 
This means that the correlated emissions occur only for 
specific directions, specific frequencies and specific polar- 
izations. There are two types of down conversions: 

1. both photons of a pair have the same polarization 
(type I), 

2. they have orthogonal polarizations (type II). 
Additionally if one has: 

1. uj ~ uj 1 then we have a frequency degenerate PDC, 

2. and if k = k' , then we have a co-linear one. 

4. Time correlations 

In this section it will be shown that the down converted 
photons are very tightly time correlated. The probability 
of a detection of a photon, of say, the horizontal polar- 
ization H, at a detector situated at point x and at the 
moment of time t, is proportional to 

p(S,t,H) ~ r ] TrQE i ~ ) (xA)E i + ) (x 1 t), (A10) 

where 77 is the coefficient which characterizes the quan- 
tum efficiency of the detection process, g is the density 
operator, which represents the field in the Heisenberg 
picture, Eh is the horizontal component of the field in 
the detector. For the above relation to be true, we also 
assume that via the aperture of the detector enter only 
photons of a specified direction of the wave vector. 

For a pure state, HAlOjl reduces to p(x, t, H) ~ 
(iplE^ E^\ip). The probability of a joint detection of 
two photons, of polarization H , at the locations x\ and 
i?2, and at the moments of time t\ and ti is proportional 
to 

p(xi 1 t 1 ;x 2l t 2 ) ~ 

(i,\E { H ] (x^E^ (x 2 ,t 2 )E { ^ ] (x 2 ,t 2 )E ( + ) (x 1 ,t 1 )\^). 

(AH) 

The evolution of the field within the crystal lasts for the 
time around — . If the detectors are very far away from 
each other, and from the crystal, then the photon field 
can be treated as free-evolving. State \4>) is the photon 
state that leaves the crystal at the moment around — : 

W - -)) = + le~ « f € H % L (t')dt'\£l). 

c in Jo 

(A12) 
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Let t = h i h = t', and \tp) = \ip(t = §)), then (jXTTT> (fi|ai(a/)at(a;)|fi) = 6ij6(w'-u>), we get 
can be written down as 



p(xi,t\x 2 ,t') ~ 

(A13) 

If we choose just two conjugate propagation directions 
(i.e. such that fulfill the phase-matching conditions), 
then approximately one has 



4 +) 

E { +\x',t' 



J duie~ 
J dwe~ 



^*/i(w)ai(w) 



where /i and / 2 are the frequency response func- 
tion which characterize the detections (or rather filter- 
detector system) . We assume that the response functions 
are such that their maxima agree with the frequencies 
determined by the phase-matching conditions. The an- 
nihilation and creation operators which were used above, 
are replaced by ai(oj) and its conjugate, which can be 
used to describe "unidirectional" excitations of the pho- 
ton field (i.e., we assume that the detectors see only the 
photons of a specified duration of propagation, this a 
good assumption if the detectors are far from the crys- 
tal, and the apertures are narrow). The index i defines 
the direction (fixed) of the wave vector. The opera- 
tors satisfy commutation relation, which are a modifi- 
cation of those given above to the current specific case 
[a; (a;), a] (a/)] = 5 10 8(lo -u'), [ai(w),a 3 -(a/)] = 0. 

Introducing an unit operator / = l^*)(^i|> where 

\bi) is a basis state into <|A13|) . we obtain 

p(xi,t\x 2 ,t') ~ 

(^\E { -\x 1 ,t)E^\x 2 ,t , )iE { +\x 1 ,t)E { +\x 2 ,t')\i>). 

(A14) 

Since E^ contains only the annihilation operators, they 
transform the two photon state |\E') into the vacuum 
state. Thus, the equation (|A14(1 can be written down 

as 



P {x x ,t\x 2 ,t') ~ (^E^E^l^iniE^E^'^), 

(A15) 

where the primed expressions pertain to the mo- 
ment of time t' and the position x 2 . Thus we 
have p(x u t\x 2 ,t') ~ \A 12 (t, t')\ 2 , where A 12 (t,t') = 

(n\E ( ^ ] {x 1: t)E^ +) {x 2l t')\^) . The state |*) can be ap- 
proximated by 

|*) = \Q)+ J duii J <kj2FQ5{w-wx-u 2 )a\(ux)a\{u) 2 )\n). 

(A16) 

Then 

A 12 (t,t') = (niJdLu'e-^'tfciLu'^iu') 

J due~ tuJt fi(ui)ai(uj) J dco 2 J duj 2 FoS(ujo — o>i — u> 2 ) 

ai,.(Lj 2 )a\{uJi)\ri) 

(A17) 

Since the creation and annihilation operators for 
different modes commute, and since one can use 



A 12 (t,t') = e ~ iu ° t ' / duFoe-^-^MtJo-^hiu;). 



(A18) 

If Fq varies very slowly, which is a good asumption in the 
case of a crystal, then we have 

p(*i,t|j?a,f)~l^ia(*,f)| 2 r . 1Q , 
~|/do; e --( t - i ')/ 2 (a;o-a;)/ 1 (a;)| 2 . [ALy) 

We see that the probability will depend on the difference 
of the detection times. 

To illustrate the above, let us assume that fi=f 2 = 
f, and that they are gaussian. Then, if one assumes 
that the central frequency of / is lo c = i^o and f(ui) — 

Ce ^ , then we have fi\u>) — f 2 \u>o — uS) = }{lo). If 
one uses these relations the probability of detection of two 
photons at the moments t and t' satisfies the following 
dependence 



p(xi,t\x 2 ,t') 
- \J dwe-^-^C 



-(t-t'f 



(A20) 

We see, that if a — > oo (that is, the detector has an 
infinitely broad frequency response) then the expression 
(|A20|I approaches the Dirac delta 6(t—t') and this means, 
that the two detectors register the two photons at the 
same moment of time. However, such detectors do not 
exist. Nevertheless, from equations <|A19() and (|A20|) . we 
see that the degree of time correlation of the detection 
of the PDC photons depends on the frequency response 
of the detectors. Thus, the photons are almost perfectly 
time correlated. 

We have shown what are the reasons for the properties 
of the PDC photons. Although the above reasoning was 
done under the assumption of a monochromatic nature 
of the pumping field, all this can be generalized to the 
non-monochromatic case, including the most interesting 
one of a pulsed pump. The distinguishing traits of this 
situation can be summarized by the following remarks. 
The emission from the crystal can appear only when the 
pulse is within the crystal. Further, the frequency ojq and 
the wave vector are not strictly defined. If the pulse is 
too short because of the relation T w x— , where T is the 
pulse width, the PDC photons are less tightly correlated 
directionally. 

The two photon state coming out of a PDC can then 
be approximated by 

|*)= fdui F (ujo) f dun f du 2 x 

A(oj - loi - uj 2 )a\(oji)al(uj 2 )\fl), 

where we have replaced the effective pump amplitude 
by the spectral decomposition of the laser pulse F(u>o). 
Since a pulse is a superposition of monochromatic waves, 
we therefore integrate over the spectrum. 
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